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HIGHLIGHT  SUMMARY 

The  capability  to  analyze  reinforced  embankments  has 
been  Incorporated  Into  the  STABL  slope  stability  analysis 
program.  The  new  program  Is  called  STABL6.  This  user's 
manual  describes  how  simplified  Bishop's  method  was  modified 
to  include  reinforcing  layers.  The  modified  program  can  an- 
alyze multiple  reinforcing  layers  and  allows  the  orientation 
of  the  reinforcing  force  to  be  specified.  Input  Instructions 
and  an  example  problem  for  the  reinforcing  layer  option  are 
gi  ven. 


PC-STABL6  WITH  REINFORCING  LAYER  OPTION 
--  USER'S  MANUAL  ~ 


INTRODUCTION 

A  slope  stability  program  PC-STABL5  developed  at  Purdue 
University  was  modified  to  give  It  the  capability  to  analyze 
embankments  constructed  on  soft  ground  reinforced  with  one 
or  more  layers  of  reinforcement.  This  modified  program  is 
called  PC-STABL6.  It  was  developed  as  part  of  a  study  of 
design  of  reinforced  embankments  (Humphrey  and  Holtz,  1986a) 
that  also  Included  a  finite  element  analysis  procedure 
(Humphrey  and  Holtz,  1986b).  PC-STABL5  is  a  general  purpose 
limiting  equilibrium  slope  stability  analysis  program  with 
solution  options  that  use  the  simplified  Janbu  method  with 
circular,  sliding  block,  or  irregular  shaped  sliding  sur- 
faces or  the  simplified  Bishop  or  Spencer  methods  with  cir- 
cular surfaces.  Program  operation  Is  described  In  Slegel 
(1975a),  Love  11,  et  at.  (1984),  and  Carpenter  (1985). 
Further  details  on  program  development  are  given  In  Slegel 
(1975b)  and  Boutrup  (1978).  The  program  Is  written  in 
FORTRAN77  and  operates  on  IBM-PC,  IBM-XT,  IBM- AT,  or  compat- 
ible computers. 

The  simplified  Bishop  analysis  routine  was  modffied  to 
Include  the  stabilizing  moment  provided  by  one  or  more 


layers  of  reinforcement.  The  user  Inputs  the  location  of 
each  layer,  the  distribution  of  available  force,  and  the 
direction  In  which  the  force  acts.  The  circular  trial  sur- 
face may  be  specified  or  PC-STABL5's  surface  generation 
routine  may  be  used. 

This  report  Is  a  user's  manual  for  analysis  of  rein- 
forced embankments  with  PC-STABL6.  Only  aspects  of  the  pro- 
gram specific  to  the  reinforcement  option  are  covered.  The 
user  Is  referred  to  the  references  noted  above  for  general 
Instructions  on  program  use.  In  the  following,  the  method 
used  to  Incorporate  the  stabilizing  effect  of  the  reinforc- 
ing layer  Is  described  and  an  outline  of  the  solution  proce- 
dure Is  given.  Appendix  A  lists  references  cited.  Appendix 
B  gives  an  example  problem.  The  data  Input  format  Is  de- 
tailed In  Appendix  C.  Error  codes  used  by  the  reinforcing 
layer  routines  are  given  in  Appendix  D.  Appendix  E  lists 
the  modifications  made  to  PC-STABL5  and  two  new  subroutines 
that  were  added.  Modifications  to  the  routine  to  output 
problem  geometry  on  a  Hewlett-Packard  7470A  two  pen  plotter 
are  listed  In  Appendix  F. 


ANALYSIS  OF  REINFORCING  LAYERS  WITH 
SIMPLIFIED  BISHOPS  METHOD 


Background 

Tensile  reinforcing  layers  provide  a  resisting  moment 
which   increases  embankment  stability.    This  effect  was 


Incorporated  Into  the  simplified  Bishop  (1955)  method  of 
analysis.  It  was  assumed  that  the  reinforcement  provides 
only  a  resisting  moment  and  does  not  alter  the  stresses  on 
the  assumed  slip  surface.  Finite  element  analyses  confirmed 
that  the  reinforcement  has  an  Insignificant  effect  on  normal 
stresses  on  the  portion  of  the  failure  surface  passing 
through  the  embankment  (Humphrey  and  Holtz,  1986a). 
Furthermore,  the  soil  strength  and  available  reinforcing 
force  are  assumed  to  be  mobilized  simultaneously.  The 
Implications  of  this  assumption  should  be  critically  exam- 
ined when  the  foundation  soils  reach  a  peak  shear  strength 
at  small  strain  followed  by  strain  softening. 

In  the  simplified  Bishop  method  a  circular  shaped  slip 
surface  is  assumed  and  divided  into  a  number  of  slices.  A 
slip  surface  and  the  forces  acting  on  a  typical  slice  are 
shown  In  Fig.  1.  The  reinforcing  force  FD  acting  at  the 
intersection  of  the  slip  surface  and  the  reinforcement  is 
also  shown.  F  has  units  of  force/unit  width.  It  provides 
a  resisting  moment  equal  to  FR  times  Its  moment  arm,  y, 
about  the  center  of  the  circle.  If  there  are  multiple  rein- 
forcing layers,  the  total  resisting  moment  is  the  sum  of  the 
resisting  moment  provided  by  each  layer.  This  additional 
resisting  moment  is  Included  In  the  equation  for  the  safety 
factor  (Ingold,  1982;  Humphrey  and  Holtz,  1986): 
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Figure  1  Simplified  Bishop  method  of  slices  Including  hori- 
zontal reinforcing  force  showing  forces  acting  on 
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where:     W    =  weight  of  the  1    slice 

Q.   =  vertical  surface  load  applied  to  the  1 
1  v 

s  1  Ice 
Q1W  =  horizontal  surface  load  applied  to  the  1 

1  n 

s  1  Ice 
d.   =  moment  arm  of  Q.. 
u.   =  pore  water  pressure  acting  on  base  of  1 

s  1  Ice 
Ax   =  width  of  1    slice 
6.   =  Inclination  of  base  of  1    slice 
c.   =  cohesion  on  base  of  I    slice 
♦.   =  friction  angle  on  base  of  1    slice 
r    =  radius  of  assumed  trial  circle 


n    =  number  of  slices 
RJ 


Fe .  =  force  in  J    reinforcing  layer 


y.   =  moment  arm  for  J    reinforcing  layer 
m   =  number  of  reinforcing  layers 
SF   =  safety  factor 

The  safety  factor  which  satisfies  Eq.  1  is  found  by  trial 
and  error.   This  equation  was  Implemented  in  STABL5. 


I mp 1 ementat i on 

To  use  the  reinforcing  option,  the  location  of  each 
layer  Is  specified  by  a  series  of  x,y  coordinates  starting 
at  the  left  end  of  the  reinforcement  and  moving  to  the  right 
as  shown  In  Fig.  2.  In  a  typical  application  the  left  end 
of  the  layer  would  be  at  the  embankment  toe  or  the  face  of 
the  slope  and  would  extend  across  the  ful  1  width  of  the 
embankment.  The  layer  may  be  horizontal  or  Inclined,  how- 
ever, the  user  must  employ  Judgment  on  the  applicability  of 
the  analysis  method  if  the  layers  are  more  than  moderately 
Inclined.  It  Is  the  user's  responsibility  to  define  a  rea- 
sonable reinforcement  geometry. 

The  available  force  and  orientation  of  the  force  Is 
specified  at  each  point  defining  the  reinforcing  layer. 
Similar  procedures  were  used  by  Duncan,  et  al.  (1985). 
Suggestions  on  choice  of  the  available  force  and  Its  orien- 
tation are  given  in  Humphrey  and  Holtz  (1986a).  The  orien- 
tation of  the  force  is  specified  by  the  Inclination  factor 
I  .  I-  varies  from  0.0  which  specifies  that  the  force  acts 
in  the  direction  of  the  reinforcement  to  1.0  which  specifies 
that  it  acts  tangent  to  the  slip  surface  as  shown  in  Fig.  3. 

The  right  most  Intersection  of  each  reinforcing  layer 
and  the  trial  surface  Is  located.  Then  the  force  <FRj)  and 
inclination  factor  (IfJ>  at  the  intersection  Is  found  by 
linear  interpolation  between  the  adjacent  specified  points 
(Fig.  2).   It  is  possible  for  the  toe  of  some  trial  circles 
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Figure  2   Location  of  reinforcement  and  distribution  of 
force  and  Inclination  factor. 
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Figure   3      Orientation   of    reinforcement    force 


to  Intersect  a  layer  a  second  time  but  this  Is  assumed  to 
have  no  effect  on  stability. 

The  moment  arm  (y)  Is  a  function  of  the  trial  circle 
radius  (r),  the  slope  of  the  reinforcement  (♦) ,  the  Inclina- 
tion factor,  the  coordinates  of  the  circle  center  (*_•  vc** 
and  the  coordinates  of  the  Intersection  (x.,  y.)  as  shown  on 
Fig.  2.  For  a  horizontal  layer  and  I  _  =  0.0,  y  =  y  -  y, 
and  for  I-  =  1.0,  y  =  r.   For  other  cases 

A  =  atan[(xc-x1 )/(yc-yf )]  (2) 

y  =  r  sln[n/2  -  A  +  *  +  If(A-*)]  (3) 

The  reinforcing  option  Is  implemented  only  for  the  sim- 
plified Bishop  method.  This  restricts  solution  options  to 
SURBIS  for  a  single  specified  trial  circle  or  CIRCL2  for 
multiple  circles  generated  using  a  random  technique.  An 
example  problem  using  the  SURBIS  option  is  shown  in  Appendix 
B.  Input  instructions  are  given  in  Appendix  C  and  error 
codes  are  listed  in  Appendix  D. 

STABL5  has  the  capability  to  analyze  many  different 
types  of  problems.  Some  combinations  of  solution  options 
may  produce  unrealistic  results,  for  example,  using  the 
reinforcing  and  tieback  options  together.  The  user  should 
exercise  good  Judgment  in  this  regard. 
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SUMMARY 

A  modification  to  a  general  purpose  slope  stability 
program  STABL5  Is  described  which  gives  It  the  capability  to 
analyze  reinforced  embankments  constructed  on  soft  ground. 
The  modification  Is  based  on  the  simplified  Bishop  method  of 
slices.  Multiple  reinforcing  layers  with  specified  distri- 
butions of  available  force  may  be  used.  The  equations  and 
procedures  used  in  the  Implementation  are  given.  An  example 
problem.  Input  Instructions,  error  codes,  and  modifications 
to  STABL5  are  given  in  the  appendices. 
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APPENDIX  B  -  EXAMPLE  PROBLEM 

This  is  a  simple  example  on  the  use  of  the  reinforcing 
layer  option  for  a  6- foot  high  embankment  with  a  2hrlv  side 
slope  as  shown  In  Fig.  B- 1  .  The  embankment  fill  Is  cohe- 
sion less  with  a  unit  weight  of  125  pcf  and  a  friction  angle 
♦  of  30° •  The  foundation  Is  soft  clay  with  a  unit  weight  of 
115  pcf  and  an  undrained  shear  strength  of  125  psf.  The 
water  table  is  assumed  to  be  very  deep.  The  embankment  is 
reinforced  with  a  single  layer  located  at  Its  base.  The 
available  force  is  specified  to  be  zero  at  the  toe  increas- 
ing to  500  lb/ft  under  the  central  portion  of  the  embankment 
as  shown  in  Fig.  B- 1  .  The  force  is  assumed  to  act  in  the 
direction  of  the  reinforcement  so  the  Inclination  factor  is 
0.0.  The  safety  factor  is  calculated  for  the  trial  circle 
shown  In  Fig.  B- 1 . 

The   following   input   file   was   prepared  using   the 

instructions  in  Slegel  (1975a)  and  Appendix  C.  The  output 

for  the  problem  Is  shown  on  the  following  pages.  A  plot  of 

problem  geometry  generated  using  PLOTSTBL.BAS  on  a  Hewlett- 
Packard  7470A  plotter  is  also  shown. 
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PROF  I L 

REINFORCED  EMBANKMENT  -  EXAMPLE 

4  3 

0.  11. 

14.  11.  2 

14.  1  1 

.  26.  17.  1 

26.  17 

.  45.  17.  1 

14.  1  1 

.  45.  11.  2 

SOIL 

2 

125.  125.  0.  30.  0.  0.  0 

115.  115.  125.  0.  0.  0.  0 

RE  INF 
i 

i 

3 

14.  1  1 

.  0.  0. 

26.  1  1 

.  500.  0. 

45.  1  1 

.  500.  0. 

SURBIS 

13 

6.00  11.00 

8.  12  8 

.88 

10.61 

7.20 

13.37 

6.03 

16.30 

5.40 

19.30 

5.33 

22.26 

5.84 

25.06 

6.90 

27.62 

8.48 

29.83 

10.51 

31  .61 

12.92 

32.90 

15.63 

33.25 

17.00 

EXECUT 
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Output  for  example 


••  PCSTABL5  "' 


by 
Purdue  University 


— Slope  Stability  Analysis- 
Simplified  Janbu,  Simplified  Bishop 
or  Spencer's  Method  of  Slices 


Run  Date: 

Time  of  Run: 

Run  By: 

1 nput  Data  F  i 1 ename : 

Output  F  i 1 ename : 

Plotted  Output  Filename: 

6/13/86 
8:30  PM 
D.  Humphrey 
exampi e. in 
exampi e.out 
exampi e.pl t 

PROBLEM  DESCRIPTION    REINFORCED  EMBANKMENT 

-  EXAMPLE 

BOUNDARY  COORDINATES 

3  Top 

4  Total 

Boundar  i  es 
Boundaries 

Boundary 
No. 

X-Left 
(ft) 

Y-Left     X-Right 
(ft)        (ft) 

Y-Rlght 
(ft) 

Soi 1  Type 
Below  Bnd 

1 

2 
3 
4 

0.00 
14.00 
26.00 
14.00 

11.00       14.00 
11.00       26.00 
17.00       45.00 
11.00       45.00 

1  1.00 
17.00 
17.00 
1  1  .00 

2 

1 
1 
2 

ISOTROPIC  SOIL  PARAMETERS 


2  Type ( s )  of  So  I  1 


Sol  1 

Type 

No. 

1 
2 


Total 

Unit  Wt 

(pcf) 

125.0 
115.0 


Saturated 

.  Unit  Wt. 

(pcf) 

125.0 
1  15.0 


Cohesion 

Intercept 

(psf ) 

0.0 
125.0 


Friction 
Angle 
(deg) 

30.0 
0.0 


Pore 

Pressure 

Param. 

0.00 
0.00 


Pressure 

Constant 

(psf) 

0.0 
0.0 


Plez. 
Surface 
No. 

0 
0 
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Output  for  example  (cont.) 
REINFORCING  LAYER(S) 

1  REINFORCING  LAYER(S)  SPECIFIED 

REINFORCING  LAYER  NO.    1 

3  POINTS  DEFINE  THIS  LAYER 

POINT      X-COORD    Y-COORD    FORCE    INCLINATION 
NO.  FACTOR 

1  14.00  11.00  0.00  0.000 

2  26.00      11.00     500.00      0.000 

3  45.00      11.00     500.00      0.000 

Trial  Failure  Surface  Specified  By  13  Coordinate  Points 


Point 

X-Surf 

Y-Surf 

No. 

(ft) 

(ft) 

1 

6.00 

1  1  .00 

2 

8.  12 

8.88 

3 

10.61 

7.20 

4 

13.37 

6.03 

5 

16.30 

5.40 

6 

19.30 

5.33 

7 

22.26 

5.84 

8 

25.06 

6.90 

9 

27.62 

8.48 

10 

29.83 

10.51 

1  1 

31  .61 

12.92 

12 

32.90 

15.63 

13 

33.25 

17.00 

Circle  Center  At  X  =    18.1  ;Y=    21.0   and  Radius,    15.7 

Factor  Of  Safety  For  The  Preceding  Specified  Surface  =   1.079 

WARNING  -  Factor  Of  Safety  Is  Calculated  By  The  Modified  Bishop 

Method.  This  Method  Is  Valid  Only  If  The  Failure  Surface 
Approximates  A  Circle. 
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Output    for    example    (cont.) 

Y                               AXIS                               FT 
0.00                 5.63               11.25               16.88              22.50              28.13 
X  0.00    + + * + + + 


5.63  + 

S 

S 
A      11.25  + 


X      16.88  + 


I      22.50  + 


28. 13  + 


33.75  + 


39.38  + 


T      45.00  + 


S 


S 


s 


s  s 
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APPENDIX  C  -  DATA  INPUT  FORMAT  FOR  REINFORCING  LAYERS 
C. 1  Input  for  Reinforcing  Layers 

COMMAND  CARD   RE  I NF     Command  Code 

DATA  CARD      Integer   Number  of  reinforcing  layers 

NOTE:   Repeat  the  following  set  of  data  cards  for 

each  reinforcing  layer. 

DATA  CARD      Integer   Number  of  points  defining  the 

reinforcing  layer 

DATA  CARD     Real      X  coordinate  of  point  on 

reinforcing  layer 

Real      Y  coordinate  of  point  on 
reinforcing  layer 

Real      Force  in  reinforcement  at  point 
(force/unit  width) 

Real      Inclination  factor;  between  0.0  and 
1.0;  =  0.0  force  acts  in  plane  of 
reinforcement;  =  1.0  force  acts 
tangent  to  failure  surface 


NOTE :    Repeat  preceding  data  card  for  each  point 
defining  the  reinforcing  layer. 
C.2  Input  for  Suppressing  of  Reactivating  Reinforcing  Layers 
COMMAND  CARD   RE  INF     Command  Code 
DATA  CARD      Integer   Number  zero  (0) 
C.3  Input  Restrictions 

1.  No  more  than  10  reinforcing  layers  can  be  specified. 

2.  A  reinforcing  layer  must  be  specified  by  at  least  2  but 
not  more  than  40  points. 

3.  The  reinforcing  force  must  be  positive. 
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The  Incl Inatfon  factor  must  be  between  0.0  and  1.0.   A 
factor  of  0.0  specifies  that  the  reinforcing  force  acts 
in  the  direction  of  the  reinforcement.   A  factor  of  1.0 
specifies  that  the  force  acts  tangent  to  the  failure 
surface. 
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APPENDIX  D  -  RE  INF  ERROR  CODES 

RF01  -  An  attempt  has  been  made  to  suppress  or  reactivate 
undefined  reinforcing  layers.  Data  must  be  defined 
by  prior  use  of  command  RE  INF  before  they  can  be 
suppressed.  Suppressed  data  can  not  be  reactivated 
If  command  PROF  I L  has  been  used  In  the  execution 
sequence  subsequent  to  use  of  RE  INF,  whether  the  data 
are  active  or  suppressed. 

RF02  -  The  number  of  reinforcing  layers  specified  exceeds 
10.  The  problem  must  either  be  redefined  so  fewer 
reinforcing  layers  are  used  or  dimensioning  of  the 
program  must  be  Increased  to  accommodate  the  problem 
as  defined. 

RF03  -  The  number  of  points  defining  a  reinforcing  layer 
exceeds  40.  The  problem  must  be  either  redefined  so 
fewer  points  are  used  or  the  dimensioning  of  the 
program  must  be  Increased  to  accommodate  the  problem 
as  defined. 

RF04  -  A  negative  coordinate  has  been  specified  for  the 
reinforcing  layer  and  point  number  Indicated.  All 
problem  geometry  must  be  located  within  the  first 
quadrant . 


RF05  -  A  negative  reinforcing  force  has  been  specified  for 
the  reinforcing  layer  and  point  number  Indicated. 
The  reinforcing  force  must  be  zero  or  positive. 

RF06  -  An  inadmissible  inclination  factor  has  been  specified 
for  the  reinforcing  layer  and  point  number  indicated. 
The  Inclination  factor  must  be  between  0.0  and  1.0. 

RF07  -  The  reinforcing  layer  and  point  number  Indicated  Is 
not  to  the  right  of  the  points  specified  prior  to  It. 
The  points  defining  the  reinforcing  layer  must  be 
specified  In  left  to  right  order. 

RF08  -  The  reinforcing  layer  Indicated  Is  specified  by  only 
one  point.  The  reinforcing  layer  must  be  specified 
by  at  least  two  points. 

RF09  -  An  attempt  has  been  made  to  use  the  reinforcing  layer 
option  with  the  simplified  Janbu  method  of  analysis. 
The  reinforcing  layer  option  can  only  be  used  with 
the  simplified  Bishop  method  of  analysis. 
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APPENDIX  E  -  MODIFICATIONS  TO  PC-STABL5 


This  appendix  describes  modifications  made  to  PC-STABL5 
which  was  originally  compatible  with  Microsoft  FORTRAN, 
Version  3.21  (MS-FORTRAN).  The  modifications  were  made  for 
two  reasons.  The  first  were  changes  to  make  the  program 
compatible  with  Ryan-McFar land.  Version  2.00  FORTRAN  (RM- 
FORTRAN).  These  were  mostly  minor  such  as  commas  In  format 
statements  and  are  listed  below.  In  addition,  MS-FORTRAN 
uses  a  special  $NOFLOATCALLS  statement  at  the  beginning  of 
each  program  unit;  these  were  made  Inactive  in  the  RM- 
FORTRAN  version  by  substituting  CNOFLOATCALLS.  The  MS- 
FORTRAN  version  writes  program  output  simultaneously  to  a 
file  (unit  6)  and  to  the  screen  (unit  •).  This  was  done 
with  two  WRITE  statements,  one  for  each  unit.  This  scheme 
is  not  compatible  with  RM-FORTRAN  since  It  treats  unit  6  and 
unit  "  as  the  same  unit.  To  eliminate  the  problem  all  WRITE 
statements  to  unit  •  were  made  Inactive  by  Inserting  *c  ••' 
In  columns  1-4.  The  second  reason  for  the  changes  was  to 
Incorporate  the  reinforcing  layer  option.  The  modifications 
are  listed  below.  They  Included  minor  modifications  to  sev- 
eral existing  subroutines  and  the  addition  of  two  new 
subroutines . 
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Modlf Icatlons  to  main  program  5TABL5 


Dana  N.  Humphrey,  Graduate  Research  Assistant,  May  1986. 


Add  after  stbl   12 

c   1986  Modifications  for  reinforcing  layer  option  (PCSTABL6) 
c 
..c 
c 
Add  after  stbl  160 

common  /blk23/lrelnf .nrelnf ,nrpts( 10) ,xrelnf( 10,40) , 

yelnf(10,40),rforce(10,40),rlnclf(10,40),rmoms 
Add  after  stbl  168 

real  load, 1 Imlt.kcoef, tncl In, length 
Integer  sltp.bn 
Replace  existing  Jrcmay85  following  stbl  172  with 

1    'LIMITS', 'ANISO','SURBIS','SPENCR','REINF'/ 
Add  after  stbl  200 
lrelnf=0 
rmoms=0. 
Replace  existing  Jrcmay85  following  stbl  324  with 

23  do  12  1=1,18 
Replace  existing  Jrcmay85  following  stbl  328  with 

1    13,30,35),! 
Replace  existing  Jrcmay85  following  stbl  350  with 
c  •••  writer, I02){keyw(l), 1  =  1, 18) 
write(6,l02)(keyw(l),l=l,18) 
Add  after  stbl  510 
35  call  relnf 
go  to  1 1 


dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 

dnhmay86 
dnhmay86 

dnhmay86 
dnhmay86 

dnhmay86 

dnhmay86 
dnhmay86 

dnhmay86 

dnhmay86 

dnhmay86 
dnhmay86 

dnhmay86 
dnhmay86 


Modifications  to  subroutine  BL0CK2 

Add  after  blk2  56 
real  limit 


dnhmay86 


Modifications  to  subroutine  EXECUT 

Add  before  exec  46 

common  /blk23/l relnf , nrelnf ,nrpts( 10) .xrelnf (10,40), 

yrelnf(10,40),rforce(10,40),rlnclf(10,40),rmoms 
Add  after  exec  62 

Ifdrelnf.eq.Dcall   re1nf2 


dnhmay86 
dnhmay86 

dnhmay86 


Modifications  to  subroutine  FACTR 

Add  after  fctr  148 

common  /blk23/lrelnf .nrelnf ,nrpts( 10), xrelnf ( 10,40), 

yrelnf( 10,40), rforce( 10,40), rlnclf( 10, 40),  rrooms 
Recpace  existing  fctr  296  with 

do  10  J*l,20 


dnhmay86 
dnhmay86 

dnhju!86 
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Replace  existing  fctr  314  with 

lf(lrelnf.eq.l)then 

bottom=sumb-rmoms/ ( rad I us'fol d) 

I f (bottom. gt.O.) then 

fnew=  sumt/bottom 

fnew=(fold+fnew)/2. 

else 

fnew=2."fold 

endlf 

else 

fnew=  sumt/sumb 

endlf 
Add  after  fctr  322 

foldx=fold 
Replace  existing  fctr  338  with 

II Ox .'Factor  Of  Safety  Calculation  Has  Gone  Through  20  Iterations' 
Replace  existing  fctr  346  and  348  with 

if(lsearc.eq.O)   wrlte(6, 108)fs,foldx 
c  •••    If(lsearc.eq.O)   writer,  I08)fs,foldx 

108  format (//l Ox, 'Factor  of  Safety  for  the  Preceding  Surface  Is  ', 

I  'Between', f6. 3,'  and',f6.3) 
Replace  existing  fctr  374  and  376  with 

wrlte(6,l08)fs,foldx 
c  •••  wrlte(M08)fs,foldx 
Add  before  fctr  386 

rmoms=0. 
Replace  existing  fctr  414  with 

110x,'•••  The  Above  Factor  Of  Safety  Is  Misleading  •••') 
Replace  existing  fctr  430  through  438  with 

II  Ox, 'The  Factor  Of  Safety  For  The  Trial  Failure  Surface  Defined'/ 
MOx.'By  The  Coordinates  Listed  Below  Is  Misleading.'/// 

I lOx, 'Failure  Surface  Defined  By', 13,'  Coordinate  Points'/// 
1 12x, 'Point' ,6x, 'X-Surf ' ,6x, 'Y-Surf ' ,/ 
M3x,'No.',8x,'(ft)',8x,'(ft)',//) 


dnhmay86 
dnhjul86 
dnhjul86 
dnhjul86 
dnhjul86 
dnhjul86 
dnhjul86 
dnhjul86 
dnhmay86 
dnhmay86 
dnhmay86 

dnhju!86 

dnhjul86 

dnhjul86 
dnhjul86 
dnhjul86 
dnhjul86 

dnhjul86 
dnhjul86 

dnhmay86 

dnhmay86 

dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 


Modifications  to  subroutine  FSPENC 


Replace  existing  fspn  726  with 

I10x,'«,«  The  Above  Factor  Of  Safety  Is  Misleading  •••') 

Replace  existing  fspn  742  through  fspn  750  with 

11  Ox, 'The  Factor  Of  Safety  For  The  Trial  Failure  Surface  Defined'/ 
IlOx.'By  The  Coordinates  Listed  Below  Is  Misleading.'/// 
IlOx, 'Failure  Surface  Defined  By', 13,'  Coordinate  Points'/// 
M2x, 'Point' ,6x, 'X-Surf ' ,6x, 'Y-Surf ' ,/ 
N3x,'No.',8x,'(ft)',8x,'(ft)',//) 

Replace  existing  fspn  780  and  fspn  782  with 
201  wrlte(6,202) 

c  •••  writer, 202) 


dnhmay86 

dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 

dnhmay86 
dnhmay86 


Hodifi cat  ions  to  subroutine  PLOT  IN 
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Add  after  plot  62 

common  /b 1 k23/ I  re  I  nf , nre I nf , nrpt s (10), xre I nf ( 1 0 , 40 ) , 

yrelnfd 0,40), rforce( 10,40) ,rtnclf( 10,40), 
Replace  existing  plot  66  with 

dimension  plttr(ll) 
Replace  existing  plot  76  with 

I      'SRF'.'TEN'.'RNF'.'NPL'.'END'/ 
Replace  existing  plot  %  with 

write(7,100)plttr(10) 
Add  after  plot  376 
c 

c   plot  reinforcing  layers,  if  applicable 
c . 

3  lf(lrelnf.eq.0)go  to  49 

write(7,100)plttr(9) 

write(7,113)nrelnf 

do  45  n=l,nreinf 

write(7,113)nrpts(n) 

x=xreinf (n, l)/scle 

y=yreinf(n,l)/scle 

write(7,115)x,y 

nn=nrpts(n) 

do  45  1=2, nn 

x=xrelnf (n, i)/scle 

y=yreinf (n,i)/scle 

write(7,115)x,y 
45  continue 
Replace  existing  plot  384  with 
49  lf(lblk.eq.0)go  to  9 


moms 


dnhmay86 
dnhmay86 

dnhmay86 

dnhmay86 

dnhmay86 

dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 

dnhmay86 


Hodifi cat  ions  to  subroutine  PLTN 

Add  after  pltn  62 

common  /blk23/lrelnf .nreinf ,nrpts( 10) ,xreinf( 10,40) , 

yrelnf(10,40),rforce(10,40),rinclf(10,40),rmoms 
Replace  existing  pltn  64  with 

dimension  pit (49, 51) ,symb(21 ),axls(9) ,scl (9) 
Add  after  pltn  66 

real  load, Incl in, length 

Integer  bn 
Add  after  pltn  68 

character  pit" I 
Replace  existing  pltn  72  with 

1    'W\'L','S','.','/'.'T','R'/ 

Add  after  pltn  366 

c 

c   position  points  defining  reinforcing  layers,  If  applicable 
c 

10  lf(lrelnf.eq.0)go  to  39 
do  36  1=1, nreinf 
nn  =  nrpts ( I ) 


dnhmay86 
dnhmay86 

dnhmay86 

dnhmay86 
dnhraay86 

dnhmay86 

dnhmay86 

dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 


Hod! ft  cat  tons  to  subroutine  PROF1L 


Modifications  to  subroutine  RANDOM 


Mod  If  I  cat  ions  to  subroutine  READER 


Mod  if teat  tons  to  subroutine  SCALER 


Modifications  to  subroutine  SLICES 


Modifications  to  subroutine  SOILWT 
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do  36  J=l ,nn 

call  postn(xrelnf(t ,j),yretnf(l,J),lx,1y) 
plt(lx-l,ly)=symb(l) 
plt(ix,iy)=symb(2l) 
36  continue 
Replace  existing  pltn  374  with 
39  do  8  l=l,nbnd 


dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 

dnhmay86 


Add  after  prof  98 

real  load, 1 imit, incl In, bn, length 
Replace  existing  prof  150  with 

1       29x, 'Purdue  University' ) 
Replace  existing  prof  172  with 

I       1  Ox, 'Output  Filename: 
Replace  existing  prof  176  with 

201  format (1  Ox, 'Plotted  Output  Filename: 


,a) 


.a) 


dnhmay86 
dnhmay86 
dnhmay86 
dnhmay86 


Replace  existing  rand  620  with 

HOx, 'Sliding  Block  Is',f6.1,//) 


dnhmay86 


Add  after  read  46 
character  m* 1 


dnhmav86 


Add  after  seal  64 

rea 1  1 oad , i  nc 1 1 n , 1 ength 
Integer  bn 


dnhmay86 
dnhmay86 


Add  before  s 1 1  c  50 
real  load 


dnhmay86 


Add  after  slwt  40 
Integer  sltp 


dnhmay86 
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Modifications  to  subroutine  TIES 

Replace  existing  ties  222  with 

105  format(/,IOx,14h",#  ERROR  -  ,a4,6h  ••••  ,5x,4hT1e  ,13,/) 


dnhmay86 


Modifications  to  subroutine  TRANS 

Add  after  tran  68 
Integer  sltp 


dnhmay86 


Modifications  to  subroutine  WEIGHT 


Relplace  existing  lines  wght  200  through  wght  208  with 

11  Ox, ••••• • 

110x,' INPUT  ERROR  -  Trial  Failure  Surface 

HOx, '••"•"•  Extends  Above  The 

II Ox, '••••••  Ground  Surface 

11  Ox, ••' 


'/ 

dnhmay86 

'/ 

dnhmay86 

'/ 

dnhmay86 

'/ 

dnhmay86 

'/) 

dnhmay86 

New  subroutine  RE  INF 

SUBROUTINE  RE  INF 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c   

c   

c 

CNOFLOATCALLS 

COMMON  /BLK01/IANGL, IBLK, IEXIT, 1CIRC, ILIMIT, IPLOT, IREAD, ISEARC, 
1  IBLK2, ISOIL, 1STR, ISURC, ISURF, ITIES, IWAT.RD.TOL, 

1  CBLK,CPLT 

COMMON  /BLK23/IREINF,NREINF,NRPTS(10),XREINF(10,40), 

YREINF(10,40),RFORCE(10,40),RINCLF(10,40),RMOMS 
DIMENSION  ERROR(8) 
CHARACTER  ERROR*4 

DATA  ERROR/'RFOl ' , 'RF02' , 'RF03' , 'RF04' , 'RF05' , 'RF06' , 'RF07' , 
'RF08'/ 


SUBROUTINE  RE  INF 


FUNCTIONS  - 

READ  THE  NUMBER  OF  REINFORCING  LAYERS. 

IF  EQUAL  TO  ZERO,  EXISTING  REINFORCING  LAYER  DATA  IS  SUPPRESSED 
OR  RECACTIVATED  IF  PREVIOUSLY  SUPPRESSED. 

IF  GREATER  THAN  ZERO,  READS,  CHECKS,  STORES,  AND  PRINTS 
REINFORCING  LAYER  DATA. 
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C   

C   READ  NUMBER  OF  REINFORCING  LAYERS 
C   

CALL  READER(DUMMY,NREIN,0) 
C   

C   CHECK  FOR  REINFORCING  LAYER  DATA  SUPPRESSION 
C   

IF(NREIN.EQ.O)  THEN 
IF(IREINF.NE.O)  THEN 
IREINF=0 
C  •••     WRITEC108) 
WRITE(6,108) 

108  FORMAT (///, 

1  OX, 'REINFORCING  LAYER  DATA  LAYER  DATA  HAS  BEEN  SUPPRESSED') 
C   

C   CHECK  IF  DATA  DEFINED  BEFORE  REACTIVATING  IT 
C   

ELSE  IF(NREINF.EQ.O)  THEN 
C  •••     WRITEClOl)ERROR(l) 
WRITE(6,101)ERROR(1) 
101     FORMAT(/,10X,14H««»«   ERROR  -  ,A4,6H   •*•*,/) 
IEXIT=I 
C   

C   REACTIVATE  SUPPRESSED  DATA 
C   

ELSE 

IREINF=1 
C  •••     WRITEC109) 
WRITE(6,109) 

109  FORMAT (///, 

1  OX, 'SUPPRESSED  REINFORCING  LAYER  DATA  HAS  BEEN  REACTIVATED') 
ENDIF 
RETURN 
ENDIF 

NREINF=NREIN 
IREINF=1 
C   _ 

C   PRINT  NUMBER  OF  REINFORCING  LAYERS 
C   

C  •••  WRITEC,103)NREINF 

WRITE(6,103)NREINF 

103  FORMATdHl, 

///.9X, 'REINFORCING  LAYER(S) ' ,//, 

13X,I2,'   REINFORCING  LAYER(S)    SPECIFIED') 
C       

C       CHECK  REINFORCING  LAYER  STORAGE  LIMIT 
C       

IF(NREINF.GT.10)THEN 
C  •••       WRITE(M01)ERROR(2) 

WRITE(6,101)ERROR(2) 

CALL  QUIT 
ENDIF 
DO  20  I=1,NREINF 
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C   READ  NUMBER  OF  POINTS  DEFINING  LAYER 
c   

CALL  READER ( DUHMY , NRPTS ( I ),0) 
C   

C   CHECK  FOR  AT  LEAST  TWO  POINTS 
c   

IF(NRPTS(I).LT.2)THEN 
C  •••     WR1TE(«,101)ERROR(8) 
WRITE(6,IOI)ERROR(8) 
IEXIT=1 
END  IF 
C   

C   CHECK  POINT  STORAGE  LIMIT 
C   

IF(NRPTS(I).GT.40)THEN 
C  •••     WRITEC,101)ERROR(3) 
WRITE(6,101)ERROR(3) 
CALL  QUIT 
END  IF 
C   

C   PRINT  REINFORCING  LAYER  DATA  HEADINGS 
C   

C  •••   WRITEC102)  I.NRPTS(I) 
WRITE(6,102)  I.NRPTS(I) 
102   FORMAT (///, 

10X, 'REINFORCING  LAYER  NO.  ',13// 

I0X.I3,'  POINTS  DEFINE  THIS  LAYER',// 

15X, 'POINT' ,5X, 'X-COORD' ,3X, 'Y-CCORD' ,3X, 'FORCE' ,3X, 

'INCLINATION'/ 

I6X, 'NO. ',36X, 'FACTOR'/) 
C   

C   READ  AND  CHECK  REINFORCING  LAYER  DATA 
C   _ 

XP=-TOL 
YP=-TOL 
DO  10  N=1,NRPTS(I) 

CALL  READER(XREINF(I,N),IDUMMY,1) 

CALL  READER(YREINF(I,N),IDUMMY,1) 

CALL  READER ( RFORCE( I, N), I  DUMMY, 1) 

CALL  READER(RINCLF(I,N),IDUMMY,1) 
c   

C   CHECK  FOR  1ST  QUADRANT  LOCATION 
C   

IF(XREINF(I,N).LE.-TOL  .OR.  YREINF(  I  ,N)  .LE.-TODTHEN 
C  •••      WRITE(M04)ERROR(4),I,N 
WRITE(6,104)ERROR(4),I,N 
104       FORMAT(/,10X,14H«»»«  ERROR  -  ,A4,6H  ••••,5X, 
'LAYER  NO.', 13, 5X, 'POINT  NO. ',13) 
IEXIT=1 
END  IF 
C   

C   CHECK  FOR  POSITIVE  REINFORCING  FORCE 
C   

IF(RFORCE(  I  ,N) . LE.-TODTHEN 
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C  •*"       WR1TEC,I04)ERROR(5),I.N 
WRITE(6,104)ERROR(5),I,N 
IEXIT=1 
END  IF 
c   

C   CHECK  LIMITS  OF  INCLINATION  FACTOR 
C   

IF(RINCLF(I,N).LE.-TOL    .OR.    RINCLF(  I  ,N)  .GT.  1  .+TODTHEN 
C  •••                WRITEC,104)ERROR(6),I,N 
WRITE(6,104)ERROR(6),I,N 
IEXIT=1 
END  IF 
c       

C   CHECK  FOR  POSITIVE  SPACING 
c   

IF(XREINF(I,N).LT.XP  .OR.  YREINF( I ,N) .LT. YP)THEN 
C  •••      WRITEC,104)ERROR(7),ItN 
WRITE(6,104)ERROR(7),1,N 
IEXIT=1 
ENOIF 

XP=XREINF(I,N) 
YP=YREINF(I,N) 

10   CONTINUE 
c   

C   PRINT  REINFORCING  LAYER  DATA 
c   

C  •••   WRITEC106)  (N,XREINF(I,N),YREINF(I,N),RFORCE(I,N),RINCLF(I,N), 
C  •••*   N=1,NRPTS(I)) 

WRITE (6, 106)  (N,XREINF(I,N),YREINF(I,N),RFORCE(I,N),RINCLF(I,N), 
N=1,NRPTS(D) 
106   FORMAT(15X,I3,3X,3F10.2,F10.3) 
20  CONTINUE 
RETURN 
END 


New  subroutine  RE  INF 2 

SUBROUTINE  REINF2 

c  

C  

C  SUBROUTINE  REINF2 

c  

c  

C 

C  FUNCTIONS  - 

C 

C      LOCATES  INTERSECTION  BETWEEN  TRIAL  SURFACE  AND  REINFORCING 

C       LAYER(S). 

C 

C      CALCULATES  RESISTING  MOMENT  PROVIDED  BY  REINFORCING  LAYER(S). 

C 

C  

c  
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C 
CNOFLOATCALLS 

COMMON  /BLK01/IANGL,IBLK,IEXIT,ICIRC,ILIMIT,IPLOT, IREAD, ISEARC, 

1  IBLK2, ISOIL, ISTR, ISURC, ISURF, ITIES, IWAT.RO.TOL, 

1  CBLK.CPLT 

COMMON  /BLK05/NSURF, SURF (100,2) 

COMMON  /BLK15/M.MB 

COMMON  /BLK22/RADIUS 

COMMON  /BLK23/IREINF,NREINF,NRPTS(10),XREINF(I0,40). 

YREINF(IO,40),RFORCE(10,40),RINCLF(IO,40),RNOMS 

DIMENSION  ERROR(l) 

CHARACTER  ERRORM 

DATA  ERROR/ 'RF09'/ 
c   

C   CHECK  THAT  MODIFIED  BISHOPS  METHOD  SPECIFIED 
c   

IF(MB.NE.1)THEN 
C  •••   WRITEClOl)ERROR(l) 
WRITE(6,101)ERROR(1) 
101   FORMAT(/,IOX,14H••••  ERROR  -  ,A4,6H  ••••,/) 
CALL  QUIT 
END  IF 

RMOMS=0. 
c   

C   CHECK  FOR  INTERSECTION  AND  CALCULATE  RESISTING 

C   MOMENT  FOR  EACH  REINFORCING  LAYER 

c   

DO  30  I=1,NREINF 
C   

C   LOCATE  INTERSECTION  BETWEEN  TRIAL  SURFACE  AND  REINFORCING  LAYER 
c   

NRPTSI=NRPTS(I) 
DO    10  J«NRPTSI,2,-1 
JUJ-1 

DO  10  K=NSURF,2,-1 
K 1 =K- 1 
CALL  INTSCT(XREINF( I ,J1 ) ,YREINF( I ,J1 ) ,XREINF( I ,J) , 

YREINF(I,J),SURF(K1,1),SURF(K1,2),SURF(K,1),SURF(K,2), 
XINT.YINT.INTS) 
IF(INTS.EQ.l)  GOTO  20 
10   CONTINUE 
c   

C   IF  NO  INTERSECTION  GO  TO  NEXT  REINFORCING  LAYER 
c   

GOTO  30 
c   

C   INTERPOLATE  REINFORCING  FORCE  AND  INCLINATION  FACTOR  AT  INTERSECTION 
c   

20   FAC*(XREINF(I,J1)-XINT)/(XREINF(I,J1)-XREINF(I,J)) 
RFORCI=RFORCE( I ,J1 )+FAC«(RFORCE( I , J)-RFORCE( I ,J1 ) ) 
RINCLI=RINCLF( I , Jl )+FAC»(RINCLF( I , J)-RINCLF( I , Jl ) ) 
C   

C   FOR  FORCE  TANGENT  TO  CIRCLE  MOMENT  ARM  EQUALS  RADIUS 
C   
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IF(RINCLI.GE.I.-TOL)THEN 

RARM=RADIUS 
ELSE 
C   

C   CALCULATE  CENTER  OF  CIRCLE 
C   

XI=SURF(1,1) 

Y1=SURF(!,2) 

X2=SURF(2,1) 

Y2=SURF(2,2) 

X3=SURF(3,1) 

Y3=SURF(3,2) 

XHALF2*(X2  +  X3)/2.0 

YHALF2=(Y2  +  Y3)/2.0 

XCNTR=((X1«»2-X2"2)«(Y3-Y2)  -  (X2"2-X3"2)VY2-Y1 )  +  (Y3-Y1) 

1  •(Y2-Y1)»(Y3-Y2))/(2.0»((XI-X2)MY3-Y2)  -  (X2-X3)» 

2  (Y2-YI))) 

YCNTR=(X2-X3)/(Y3-Y2)«(XCNTR  -  XHALF2)  +  YHALF2 
C   

C   CALCULATE  MOMENT  ARM  FOR  FORCE  NOT  TANGENT  TO  CIRCLE 
C   

A1=ATAN((YREINF(I,J)-YREINF(I,J1))/(XREINF(I,J)-XREINF(I,J1))) 
A2=ATAN((XINT-XCNTR)/(YCNTR-YINT)) 
A3=A1+RINCLI#(A2-A1) 
A4= 1. 570796327- A2+A3 
RARM=RADIUS»SIN(A4) 
END  IF 
C   

C   CALCULATE  AND  SUM  RESISTING  MOMENT 
C   

RMOM=RFORCI#RARM 

RMOMS=RMOMS+RMOM 
30  CONTINUE 
RETURN 
END 
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APPENDIX    F    -    LIST   OF    ADDITIONS    AND    MODIFICATIONS    TO 

PLOTSTBL.BAS 


The  modifications  which  were  made  to  PLOTSTBL.BAS  to 
give  It  the  capability  to  plot  the  location  of  reinforcing 
layers     are      listed     below.  An      'A'       In     the      first      column 

Indicates  that  the  line  should  be  added  to  the  existing 
program  and  '  R'  indicates  that  the  line  should  replace  the 
exl st 1 ng    1 i  ne. 

Modifications  to  PLOTSTBL.BAS 

R  530   DIM  PLTI(ll) 

R   730   FOR    1=1    TO    11 

R    1090  FOR  J=l    TO   11 

R    1110  ON  J  GOTO    1140,1360,1570,1750,1940,2152,2310,2450,2131,2600,2700 

A  2131    REM     

A  2132  REM   PLOT  REINFORCING  LAYERS,  IF  APPLICABLE 

A  2133  REM  

A  2134  IF  PROMPT$="y"  THEN  GOSUB  2900 

A  2135  INPUT  #2,NREINFl 

A  2136  FOR  1=1  TO  NREINFX 

A  2137  INPUT  #2,NNl 

A  2138  GOSUB  2830 

A  2139  PRINT  #1,  "PU.PA" ;X#, Y#; 

R  2140  PRINT  #1,  "SI. 15, .2" 

A  2141  PRINT  #1.  "CP-.5,0;DI;LBR"+CHR$(3) 

A  2142  PRINT  #1,  "PU,PA";X#, Y#; 

A  2143  FOR  J=2  TO  NN1 

A  2144  GOSUB  2830 

A  2145  GOSUB  2770 

A  2146  PRINT  #1,  "CP-.5,0;DI ;LBR"+CHR$(3) 

A  2147  PRINT  #1,  "PU,PA";X#,Y#; 

A  2148  NEXT  J 

A  2149  NEXT  I 

R  2150  GOSUB  2770 

A  2151  GOTO  1130 

A  2152  REM  

A  2153  REM   PLOT  SEARCH  BOXES  FOR  BLOCK,  IF  APPLICABLE 

R  2950  FOR  JJ=1  TO  9 

R  2970  ON  J  GOSUB  3000,3030,3060,3090,3120,3150,3180.3210,3241 

A  3241  INPUT  "SELECT  PEN  FOR  PLOTTING  REINFORCING  LAYERS  (I  OR  2)  ";XX 

A  3242  GOSUB  3250 

A  3243  RETURN 

A  3401  DATA  RNF 


